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ABSTRACT 

The Hubble morphological sequence from early to late galaxies corresponds to an increasing rate of 
specific star formation. The Hubble sequence also follows a banana-shaped correlation between 24 and 
70 micron luminosities, both normalized with the K-band luminosity. We show that this correlation 
is significantly tightened if galaxies with central AGN emission are removed, but the cosmic scatter 
of elliptical galaxies in both 24 and 70 micron luminosities remains significant along the correlation. 
We find that the 24 micron variation among ellipticals correlates with stellar metallicity, reflecting 
emission from hot dust in winds from asymptotic giant branch stars of varying metallicity. Infrared 
surface brightness variations in elliptical galaxies indicate that the K - 24 color profile is U-shaped 
for reasons that are unclear. In some elliptical galaxies cold interstellar dust emitting at 70 and 160 
microns may arise from recent gas-rich mergers. However, we argue that most of the large range of 
70 micron luminosity in elliptical galaxies is due to dust transported from galactic cores by feedback 
events in (currently IR-quiet) active galactic nuclei. Cooler dusty gas naturally accumulates in the 
cores of elliptical galaxies due to dust-cooled local stellar mass loss and may accrete onto the central 
black hole, releasing energy. AGN-heated gas can transport dust in cores 5-10 kpc out into the hot gas 
atmospheres where it radiates extended 70 micron emission but is eventually destroyed by sputtering. 
This, and some modest star formation, defines a cycle of dust creation and destruction. Elliptical 
galaxies evidently undergo large transient excursions in the banana plot in times comparable to the 
sputtering time or AGN duty cycle, 10 Myrs. Normally regarded as passive, elliptical galaxies are the 
most active galaxies in the IR color-color correlation. 



1. INTRODUCTION 

In a series of recent papers we describe observations 
and interpretations of far infrared emission from dust in 
early-type galaxies, emphasizing emission from dust that 
is in thermal contact with hot gas virialized in galactic 
or group potential wells (Temi, Brighenti & Mathews 
2005; 2007a,b; 2008; 2009a,b = TBM05; TBM07,a,b; 
TBM08; TBM09a,b). We propose that interstellar dust 
is ejected and stripped from asymptotic giant branch 
(AGB) stars as they orbit through the hot interstellar 
gas (T w 10*^ - lO^K). Interstehar dust is heated both by 
diffuse galactic starlight and by impacts of thermal elec- 
trons. The loss of thermal energy by electron-grain im- 
pacts can cool gas near black holes in the galactic cores, 
possibly stimulating energetic feedback events (Mathews 
& Brighenti 2003). Dust both influences and is influenced 
by AGN events. In addition, mid infrared circumstellar 
emission observed in elliptical galaxies is emitted by hot- 
ter dust associated with outflows from mass-losing AGB 
stars. 

When mid (24/zm) and far (70/zm) infrared luminosi- 
ties from galaxies of all types are compared and normal- 
ized by the if-band luminosity, L24/LK and Liq/Lk, 
galaxies of all Hubble types occupy a remarkably well- 
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defined banana-shaped correlation. Remarkably, the 
galactic Hubble sequence is ordered along the correlation 
from early to late morphologies (TBM09b). For spiral 
and irregular galaxies this color-color correlation results 
from increasing specific star formation toward later type 
galaxies. Evidently, dust heated to a range of tempera- 
tures by young stars emits mid and far infrared radiation 
in nearly the same proportion, independently of the spe- 
cific star formation rate. 

By contrast, early-type elliptical galaxies, in which star 
formation is small or absent, are contiguous with the ear- 
liest spiral galaxies in the color-color plot and occupy an 
extended region of decreasing \og{L'jQ/ Lk) with nearly 
constant \og{L2A/ Lk) ~ 30.2. In general, mid and far 
infrared emission from dust grains in early-type galaxies 
does not indicate the star formation rate. 

Here we discuss an improved, sharper version of the ba- 
nana color-color plot for normal galaxies in which galax- 
ies with concentrated mid infrared emission (probably 
mostly AGNs) are removed. In addition, we propose pos- 
sible explanations for the cosmic scatter among elliptical 
galaxies in the infrared color-color plot. 

2. THE BANANA PLOT 

The top panel of Figure 1, taken from TMB09b, shows 
the banana-shaped {L2i/ Lk)-{Ljo/Lk) plot for nearby 
galaxies observed in the Spitzer Infrared Nearby Galaxies 
Survey (SINGS; Kennicutt et al. 2003). In Figure 1, and 
subsequently, L24 and L70 are in cgs units and Lk is in 
solar units. The infrared Spitzer band luminosities at 24 
and 70/im are defined by L\ = XFxAttD^ where D is the 
distance. The SINGS sample was selected to represent 
nearby normal galaxies of all morphological types. The 



2 



Mathews et al. 



banana-shaped plot in the top panel of Figure 1 provides 
an infrared main sequence for the Hubble elassification. 
Colors in Figure 1 (top) designate five morphological 
types using the de Vaucouleurs T parameter from Hy- 
perLeda (Paturel et al. 2003). The morphological type 
varies rather smoothly along the banana from early types 
at the lower left to the latest galaxies at the upper right: 
E — >• Sa,Sab Sb,Scd — )• Irr. Proceeding upward in Fig- 
ure 1, the luminosity ratio L2i/LK oc {L^o/LKy "^ is a 
surrogate for an increasing specific star formation rate 
(SFR) from early to late types (Calzetti et al. 2007). SO 
galaxies have SFRs that vary from essentially zero, simi- 
lar to red and dead E galaxies, to rates typical of irregular 
galaxies having the highest specific SFRs (TBM09a), i.e. 
SOs can be found all along the banana. 

The galaxy sample discussed by TBM09b, selected 
from the Spitzer archive, is rich in early type galaxies, 
complementing the emphasis on later types in the SINGS 
data. The central panel of Figure 1, which combines 
these two data sets, shows that the lower region of the 
banana occupied by E and E-SO galaxies is quite exten- 
sive and possible reasons for this extension are discussed 
below. 

Recently we discovered that many of the outlying 
galaxies in the TBM09b sample displaced above the ba- 
nana (indicated with green squares in Figure 1) exhibit 
strongly concentrated 24/xm emission. This concentra- 
tion is apparent from the visibility of Airy rings in 24/im 
images, including additional angular artifacts character- 
istic of MIPS diffraction at this wavelength (see also 
Calzetti et al. 2010). While it is possible that 24/xm 
diffraction in some galaxies may be due to centrally con- 
centrated starbursts, most of the galaxies with Airy rings 
have been previously identified as having additional in- 
dependent evidence of AGN activity. Moreover, in the 
central panel of Figure 1 it is seen that galaxies hav- 
ing Airy rings (green squares) lie above the banana cor- 
relation almost everywhere along the banana. If con- 
centrated starbursts are responsible for Airy rings, the 
star formation trajectory in Figure 1 of dust heated 
in unresolved starbursts would need to be steeper than 
the locus of star formation along the upper banana, 
S \og{L24 / L k) / S \og{Lio / L k) ~ 1-2, which seems un- 
likely. For these reasons, we argue that harder AGN radi- 
ation, not concentrated starbursts, is the most likely ex- 
planation for the 24/im Airy rings. All early-type galax- 
ies with Airy diffraction probably have low-level active 
nuclei (Tang et al. 2009) 

The top row of three panels in Figure 2 illustrates 24/xm 
images of two galaxies with visible Airy rings, NGC 1386 
and NGC 0315, and one, NGC 4472, which has no de- 
tectable Airy diffraction. Galaxies in Figure 1 with sig- 
nificant visible Airy rings are listed in Table 1. Cen- 
trally normalized surface brightness profiles are shown 
in the lower panel of Figure 2 together with the Spitzer 
point response function at 24/im, which decreases to half 
its central value at about 3.4 arcseconds. When (AGN) 
galaxies with strong 24//m Airy rings are removed, as 
in the bottom panel of Figure 1, the infrared color-color 
banana plot becomes much tighter. 

2.1. Bottom, of the Banana 

Figure 3 shows an expanded view of galaxies that oc- 
cupy the bottom, rather flat part of the banana, most or 



all of which are E or E-SO galaxies (T < —2.7). Galaxies 
with strong central point source emission at 24/im have 
been removed. Of particular interest is the large extent 
of the region occupied by early type galaxies, having an 
intrinsic scatter that exceeds the errors of observation 
(e.g. as shown in the vertical error bars in Fig. 4). 
If these elliptical galaxies are dominated by old stellar 
populations, one might imagine that they would occupy 
a much smaller region in the banana plot. We now dis- 
cuss in turn the physical origin of the vertical {L24/LK) 
and horizontal {L'jq/Lk) scatter of early- type galaxies in 
Figure 3. The outlying galaxy with the lowest L^aILk 
in Figures 1 and 3 is NGC 4406. 

3. INTRINSIC SCATTER OF L-za/Lk FOR E GALAXIES: 
STELLAR METALLICITY 

Low mass, oxygen-rich asymptotic giant branch (AGE) 
stars are thought to dominate the optical luminosity and 
dust production in elliptical galaxies (e.g. Habing 1996; 
Athey et al. 2002). Consequently, it is natural to ex- 
pect that AGB stars having higher metallicity also have 
more hot circumstellar dust and more emission at 24/um, 
possibly accounting for the vertical scatter in Figure 3. 
To test this idea, we use the sample of elliptical galax- 
ies with metallicities [Z/H] determined by Trager et al. 
(2000a, b), most of which have been observed at 24/im 
and which are identified with red circles in Figure 3. Rel- 
evant properties of Trager's sample are listed in Table 2. 

Figure 4 shows that L24/LK does indeed correlate with 
stellar metallicity [Z/H] = \og[{Z / H)gaiaxv / {Z / H)q] is 
normalized to the solar abundance: 

\og{L2i/LK) = 29.99 ± 0.31 + (0.52 ± Q.IQ)[Z/H] (1) 



log(L24/ix) = 30.12 ± 0.42 + (0.76 ± 0.21)[Z/iJ] (2) 

Values of [Z/H] in equations (1) and (2) refer respec- 
tively to apertures of size i?e/8 and -Re/2 where Re is 
the effective radius in the B-band (Trager et al. 2000b) . 
The correlations in Figure 4 indicate that broadband 
observations at rest-frame 24/im could be used to de- 
termine the stellar metallicity of distant, old-population 
early- type galaxies at higher redshift, provided there are 
no strong AGN contributions. The correlation between 
stellar age and metallicity among elliptical galaxies dis- 
cussed by Terlevich & Forbes (2002), is not apparent in 
the galaxies plotted in Figure 4. 

3.1. Metallicity Gradients and 24lJ,m Surface Brightness 

Proceeding further with this idea, can the (negative) 
stellar metallicity gradients in elliptical galaxies be mea- 
sured from radial variations in the K — 24/im color? 
An initial response to this question appears in Temi, 
Brighenti, Mathews (2008 = TBM08), and we follow the 
same data reduction procedure here. In TBM08 we dis- 
cussed the mean {K — 24) color profile averaged over 
a sample of only 19 galaxies, 10 of which are in com- 
mon with the Trager et al. (2000b) sample. In that 
paper we found that {K — 24) (R) decreased linearly with 
[Z/H] to about R « 1.2Re{K), then flattened, and at 
R > 1.8Re{K) began to rise again. Re{K) is the K-band 
half-light radius. Now with the complete Trager sample, 
about twice as large, we find a similar mean {K — 24) 
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differential surface brightness profile shown in the upper 
panel of Figure 5. 

The flattening of {K - 24) at 0.8 <, R/Re{K) 1.6 
and subsequent rise at larger R/Re{K) in Figure 5 does 
not indicate that these elliptical galaxies have bright ha- 
los of 24/xm emission. Figure 1 of TBM08 shows that the 
surface brightness profiles of typical elliptical galaxies at 
all Spitzer infrared wavelengths (including 24y^m) deviate 
rather little from the de Vaucouleurs profiles in optical 
light. But systematic radial color variations do exist and 
{K—2A) is the most unusual. In the lower panel of Figure 
5 we convert {K — 24) into \Z/H] using the linear rela- 
tion shown in Figure 4 with [Z / H] evaluated in the i?e/2 
aperture. We assume that the linear behavior in equation 
2 can be extrapolated to larger [Z/H]. Values of [Z/H] 
at i?e/2 are preferred to match the global luminosities 
L24 and Lk in Figure 4 and because R > Re{K) for all 
Z/H values in Figure 5. When [Z/H] values observed 
at Re /I are used, the U-shape of the radial variation in 
the lower panel of Figure 5 is largely unchanged but it is 
normalized to somewhat larger values of [Z/H]. 

Optical observations of [Z/H] profiles using the Lick 
indices at comparable R/Re{K) invariably indicate 
monotonically decreasing stellar metallicity with galac- 
tic radius and this trend might also be expected in the 
hot dust emission. But how can the flattening and pos- 
sible rise in the [Z/H] profile in Figure 5 be explained? 
If the outer regions of elliptical galaxies are formed by 
the accretion and disruption of smaller galaxies, then we 
would expect lower metallicity at large R/Re{K) typical 
of smaller galaxies, apparently opposite to the trend in 
Figure 5. Excess {K — 24) at large R/Re{K) could in 
principle be due to dust acquired by accretion in distant 
regions of the galactic atmosphere where the gas density 
(and possibly temperature) is much lower, increasing the 
grain sputtering time, preserving a higher dust density. 
However, in normal elliptical galaxies emission at 24/xm 
is dominated by hot dust in circumstellar regions around 
mass-losing AGB stars, not by emission from much colder 
interstellar dust grains, (see for example our computed 
models in Fig. 8). 

Does dust in low metallicity AGB winds have higher 
24/xm luminosities? Computed evolutionary tracks of 
AGB stars exhibit increasing effective temperature with 
decreasing abundance (Cristallo et al. 2009). If 24/xm- 
emitting grains are heated by thermal impacts, this in- 
crease in effective temperature may be sufficient to ac- 
count for the enhanced circumstellar 24/im emission at 
large R/Re{K) seen in Figure 5. Another possible ex- 
planation is that the outflow velocity in AGB stars is 
reduced (or rendered intermittent) by lowered radia- 
tion pressure acting on fewer or smaller grains in low- 
metallicity AGB stars, increasing the gas density in the 
AGB outflow, possibly allowing a smaller mass fraction 
of dust grains to emit disproportionally more in the mid- 
infrared. If low metallicity is the source of the apparent 
[Z/H] upturn in Figure 5, it must occur at metallicities 
lower than those plotted in Figure 4 where no low-[Z/H] 
upturn is visible. 

In any case, the decreasing {K — 24) at small radii, 
R/Re{K) < 0.8, in the upper panel of Figure 5 is 
consistent with the expected decrease in stellar abun- 
dance observed optically in this region. The stellar 
metallicity gradient in Figure 5 at R/Re{K) ^ 0.8, 



A[Z/ H]/ A\ogR « —0.23, is in good agreement with pre- 
viously observed metallicity gradients in elliptical galax- 
ies (CaroUo et al. 1993; Davies et al. 1993; Kobayashi 
& Arimoto 1999; Mehlert et al. 2003; Kobayashi 2004; 
Koleva et al. 2011). 

4. INTRINSIC SCATTER OF Ln/LK- COLD 
INTERSTELLAR DUST 

In the upper part of the banana, where logL^Q/Lx ^ 
31.1 in Figure 1, 24 and 70 micron emission from nor- 
mal galaxies are dominated by star formation. For 
log Ljo/Lk ^ 31.1, the 70 micron emission from early- 
type, mostly E galaxies is a measure of the vary- 
ing amounts of cold dust in the hot interstellar gas 
(TBM07a,b). 

To discuss horizontal variations and excursions in 
L-jo/Lji along the bottom of the banana, we return to 
the complete galaxy sample shown in the lower panel of 
Figure 1 and consider only elliptical galaxies. Figure 6 
shows the huge range in far infrared luminosity L70 for 
elliptical galaxies (T < —3.5) having similar Lk- Galax- 
ies with the largest Lyo have either undergone d\ist-rich 
mergers (such as NGC 5018) or are central galaxies in 
groups, although some galaxies with the lowest L70, such 
as NGC 1399, are also group-centered. L70 luminosities 
of NGC 1399 and 4472, both near the Spitzer detection 
limit, are in excellent agreement with our steady state 
model for the minimum far infrared emission from inter- 
stellar dust ejected from AGB stars into the hot ISM of 
normal elliptical galaxies as described by TBM07a,b. In 
this steady state calculation, which we now briefly re- 
view, stellar dust is created locally at the same rate that 
it is consumed by sputtering. 

Approximate stellar ages and main sequence turnoff 
rates in elliptical galaxies can be obtained from Balmer 
line absorption in optical spectra and photometry. Prom 
this follows (1) the mean current stellar mass loss rate 
from old stellar populations, dM^/dt sa l.b{M^/W^MQ) 
Mq yr~^, and (2) the associated interstellar dust pro- 
duction rate by AGB stars, assuming a metallicity- 
dependent (dust mass) /(gas mass) ratio normalized to 
~ 0.01 at solar abundance. Warm, dusty gas is deposited 
into the hot (T ~ 10^ — lO^K) interstellar gas along the 
ram-pressure stripped wakes of orbiting AGB stars and 
is subsequently maintained at T lO^K by UV from 
post- AGB stars. The stripped warm gas must thermally 
merge into the hot interstellar gas (T ~ 10^ — 10^ K) 
in less than about lO'''"'' yrs, otherwise optical Ha emis- 
sion line luminosities would exceed typical observations 
(Mathews 1990; Mathews & Brighenti 2003). 

Consequently, after ~ 10^-^ yrs dust grains can be ex- 
pected to come into direct thermal contact with the hot 
gas and begin to be eroded (sputtered) by collisions with 
H+ and He++ ions. A typical grain sputtering lifetime is 
~ 10^ yrs at galactic radius r ^ 5 — 10 kpc. Interstellar 
dust grains are heated in comparable amounts by ab- 
sorption of local starlight and by inelastic collisions with 
thermal electrons in the hot gas. Each colliding electron 
deposits ~ (3/2)fcT into the grain. 

We assume that dust ejected from AGB stars, has a 
typical initial grain size distribution No{r,a) oc a^^'^ 
normalized to the local stellar metallicity at galactic ra- 
dius r. In a steady state, when grain heating-cooling and 
formation-destruction balance at every galactic radius r, 
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we can determine the temperature T(r, a) of (amorphous 
sihcate) grains of radius a. Integrating over the grain size 
distribution (ahered by sputtering) and the projected 
surface of the galaxy, we can calculate (TBM07a,b) the 
total specific luminosity or flux XFx in any field of view 
aperture. By this means the infrared SED can be esti- 
mated for any elliptical galaxy based on its luminosity, 
distance and the temperature and density profiles of hot 
interstellar gas observed in the X-ray. 

We find that the far infrared luminosity observed in 
certain galaxies - such as NGC 1399 and NGC 4472 - 
are very close to our steady state predictions (TBM07a). 
These galaxies are also near the lower limit of ellipti- 
cal galaxies that can be detected with Spitzer MIPS at 
70/xm. For these galaxies, the local hot interstellar gas is 
assumed to consume dust by sputtering at the same rate 
that it is produced in local AGB outflows. This model 
can reproduce 70/xm and 160/zm luminosities for ellipti- 
cals with rather low Ljq, but trolly interstellar dust is 
too cold to emit observable 24/im emission. Evidently, 
L24 is produced exclusively by much hotter dust in AGB 
outflows, heated only by the nearby parent AGB stars, 
and this emission is not considered in our dust evolution 
model. 

If interstellar L^o emission from all elliptical galaxies 
agreed with our steady state predictions, they would oc- 
cupy a disconnected island near log L^q/Lk ~ 29.0 ±0.5. 
And observed AGB-heated dust would still dominate 
24/Ltm emission from ellipticals with some scatter for vari- 
able metallicity, logL^i/ Lk ~ 30.10 ± 0.15. Our steady 
state model predicts that the banana plot would appear 
only as this island for ellipticals and a separate powerlaw 
sequence of star-forming galaxies of later types. 

However, and to our initial surprise, many optically 
normal E galaxies - e.g. NGC 4636 and 5044 - have 
that are 10-30 times brighter than our steady state 
prediction. This infrared excess indicates the presence 
of excess cold interstellar dust (of mass ~ 10^ Mq) that 
cannot have come from local mass-losing stars. In addi- 
tion, elliptical galaxies with excess Ljq often have spa- 
tially extended far infrared emission, extending beyond 
the Spitzer 70/im point response fmiction (TBM07a). 
Galaxies with observationally confirmed extended dust 
are indicated in Figure 6 with large red circles. The far 
infrared color Ligo/i^To of observed excess dust emission 
conflrms that this excess dust is located further out in 
the E galaxy interstellar atmosphere, at r ~ 5 — 10 kpc, 
where grain heating by starlight and electron-grain col- 
lisions are diminished (TBM07b). Because of the short 
dust sputtering lifetime, ~ 10'' yrs, the excess dust emis- 
sion at these large radii must be transient. 

Transient, spatially extended dust can arise (1) from 
recent mergers or (2) by internally produced dust that 
is hydrodynamically transported outward from galactic 
cores with AGN-heated gas. Some elliptical galaxies with 
the highest X70 have acquired dust by a recent merger 
and sometimes a post-merger dust distribution is clearly 
visible. In other cases it is difficult to determine if dust is 
supplied externally for example by minor mergers since 
these may be difficult to detect without a detailed anal- 
ysis, as argued by Crockett et al. (2011). Nevertheless, 
our impression from images of many elliptical galaxies 
with high is that they are not surrounded by small, 



gas-rich galaxies that are about to merge. If the dust 
sputtering lifetime ^ 10^ yrs is relevant, such mergers 
would have to be quite frequent, requiring multiple simul- 
taneous ongoing small mergers to explain for example the 
extended 70^111 image of NGC 5044 (TMB07b). At a few 
kpc radius, the orbital time of merging galaxies is likely 
to be longer than the sputtering time of dust stripped 
from them. In addition, gas-rich merging galaxies must 
release widespread dust into the hot gas atmosphere that 
emits 70/Ltm, but star-forming regions in gas-rich merging 
galaxies must not emit detectable 24/Ltm emission since 
emission at this wavelength correlates with the stellar 
metallicity of the host galaxy (Fig. 4). Moreover, dwarf 
galaxies observed near central galaxies in groups/clusters 
tend to be free of interstellar gas, due in part to ram strip- 
ping in the distant past. For these reasons we consider 
here the hypothesis that internally produced interstel- 
lar dust can explain many elliptical galaxies with excess 
dust often extending far beyond the effective radius of 
the central galaxy. In TBM07b we described how dusty, 
hot interstellar gas can be buoyantly or hydrodynami- 
cally transported outward as a result of AGN heating 
events (feedback) in galaxy cores. 

For internally produced dust to explain the extended 
excess infrared emission in elliptical galaxies, it is es- 
sential to have a reservoir of dust near central black 
holes. This source of dust is conveniently provided by 
optically thick clouds or disks of dense dusty gas ob- 
served in projection against starlight in galactic cores 
(e.g. van Dokkum & Franx 1995; Lauer et al. 2005). 
The common incidence of dusty disks or chaotically dis- 
tributed dusty clouds within the central kpc of many 
ellipticals indicates that the steady state dust creation- 
destruction assumption we use to compute extended far- 
IR dust emission breaks down in dense elliptical cores. 
Apparently more dust is created in elliptical cores than 
can be destroyed by local sputtering. Dust in stellar 
ejecta can radiatively cool ambient hot gas faster than 
the sputtering or freefall time, as discussed by Math- 
ews & Brighenti (2003). When the central black hole 
accretes some nearby dust-cooled gas and releases feed- 
back energy, we propose that some of the central dusty 
gas is heated and flows buoyantly out to 5-10 kpc from 
the center where it thermally merges with the hot inter- 
stellar gas before the dust is sputtered away. We identify 
these AGN-driven outflows as the source of excess ex- 
tended far infrared emission observed in many elliptical 
galaxies shown with red circles in Figure 6. 

Figure 7 shows elliptical galaxies (T < —3.5) in the 
bottom of the banana plot in a region similar to Figure 
3. Galaxies with red dots in Figure 7 are chosen because 
of their large L^q/Lk for given L2i/LK-, lying along the 
rightmost boundary of the lower banana. These same 
galaxies are also plotted with red dots in Figure 6 where 
it is seen that they all have significant excess 70/im emis- 
sion. Galaxies with red circles in both Figure 6 and 7 are 
those in which the excess far-IR emission is known to be 
spatially extended in 70/xm Spitzer images (TBM07a,b). 

Small open green circles in Figures 6 and 7 show the 
approximate computed locations of two galaxies - NGC 
5044 and NGC 4636 using our steady-state dust evo- 
lution model in which dust ejected from galactic stars 
is destroyed locally by sputtering. The green arrows il- 
lustrate the large additional (excess) contribution to L^q 
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which we interpret as dust in gas heated by AGN feed- 
back events and then transported outward from galac- 
tic cores. Without this extended excess dust these two 
galaxies may not have been detected with Spitzer. Simi- 
lar displacements probably occurred recently in the other 
red elliptical galaxies in the banana boundary region of 
Figure 7. 

If our interpretation is correct, group-centered galax- 
ies, like NGC 5044 and NGC 4636, move back and forth 
along the bottom of the banana in Figure 1 on timescales 
determined by the grain sputtering lifetime and the inter- 
val between AGN outbursts, which may be roughly com- 
parable, ~ 10^ yrs. Elliptical galaxies that are not known 
to be group-centered can also participate in similar large, 
transient horizontal excursions in the banana plot, pro- 
vided their virialized hot gas atmospheres are hot enough 
(T > 4 X 10^ K) to sputter and are sufRciently extended. 
However, there appears to be no correlation between ex- 
cess 70/im emission and X-ray luminosity (TBM07a). 

One potential difficulty with our hypothesis is that no 
galaxy is observed to lie appreciably further to the right 
than the group of red dot galaxies in Figure 7. Compar- 
ing Figure 7 and the bottom panel of Figure 1, it is seen 
that the slope of the rightmost envelope of red galaxies in 
Figure 7 lies just along the star-formation track defined 
by galaxies of later type {{L24/LK) oc {Liq/ Lk)^''^). 
Somehow, the rightward excursions caused by the out- 
ward movement of excess cold dust due to feedback 
events must be limited so as not to exceed the star- 
formation locus. 

To understand this, we conclude that central star for- 
mation must eventually occur when the mass of dusty 
gas accumulates beyond a certain point in the cores of 
elliptical galaxies. Star formation causes both L70 and 
L24 to increase as elliptical galaxies rise along the SINGS 
star-forming locus. Star formation and AGN-driven dust 
transport into the hot gas atmospheres are not mutually 
exclusive. If central star formation commences, the hori- 
zontal green arrows in Figure 7 should bend at the right 
to move up along the star-formation trajectory with in- 
creasing L2a/ Lk- It is likely that low- level star forma- 
tion occurs in all elliptical galaxy cores (Ford & Bregman 
2012) in dense dusty clouds formed by dust-assisted cool- 
ing of gas ejected from AGE stars (Mathews & Brighenti 
2003). 

A quasi-steady-state cyclic dust evolution that involves 
reservoirs of increasing dust mass in galaxy cores may 
explain the large variation in L^q/Lk shown in Figures 
6 and 7. Cold gas clouds containing optically visible 
dust slowly accumulate in galaxy cores as gas from stel- 
lar mass loss is cooled by dust. Occasionally these dusty 
clouds may be heated by AGN events to ^ 10^ K and 
flow out into the surrounding hot atmosphere emitting 
70mum radiation. Ultimately this extended dust is likely 
to be sputtered into the gas phase. Such a hypothetical 
dust evolution could create and destroy dust at approxi- 
mately balancing rates. Star formation may or may not 
occur during each feedback episode, and it would also 
serve to consume some of the accumulated cold gas and 
dust. 

The evolution of optically thick dust clouds in the cores 
of elliptical galaxies must also participate in the quasi- 
cyclic production of extended cold dust. If cold dusty gas 



within the central ~kpc forms by cooled stellar ejecta as 
we propose, the default endstate is most likely a small ro- 
tating dusty disk. The spin axis of the disk will be aligned 
with the kinematical minor axis of nearby stars. While 
such central disks are often observed (e.g. van Dokkum 
& Franx 1995; Lauer et al. 2005), it is more common to 
find a multitude of patchy, chaotically distributed dust 
clouds or no dust absorption at all. Because of the small 
solid angle presented to the central black hole, thin, cold 
dusty disks may be difficult to disrupt if AGN energy is 
released only at the black hole. However, disks may be 
more easily disrupted by symmetric depositions of AGN 
energy in bipolar regions somewhat removed from the 
black hole and misaligned from the disk axis. In this case 
shock energy reaches opposite sides of the disk at differ- 
ent times, causing mechanical disk disruption. Even if 
strong shock waves arrive simultaneously on both sides, 
if the shocks are strong enough, much or all of the disk 
could be sufficiently heated to merge with the hot gas. 
While these important details remain unexplored, in at 
least one case we observe both a relaxed dusty disk with 
nearby dense clouds presumably from a disrupted pre- 
vious disk generation. In NGC 5044, where extended 
dust is resolved at 70/im (TMB07b), the central kpc con- 
tains many chaotically distributed optically thick clouds 
with orbital lifetimes ~ 10^ yrs compatible to the hydro- 
dynamic and sputtering lifetimes of the extended dust 
(TBM07b). In many cases remnants of disrupted central 
dusty gas may be too small or too widely dispersed by 
AGN outbursts to be visible as optically thick clouds. 

5. UN-NORMALIZED BANANA PLOT 

Figure 8 shows a less compact version of the banana 
plot in which L24 and L^q are plotted without normaliza- 
tion with Lk- sings galaxies (crosses in Fig. 8) define 
a star-forming sequence L24 oc L^q^''^^ in which ^24/^70 
is approximately constant as ^24 varies by five orders of 
magnitude. Many of the Spitzer-archived elliptical galax- 
ies in Figure 8 that occupy the SINGS star-forming re- 
gion above logL24 ~ 41.8 are pecufiarQ SINGS galaxies 
in Figure 8 having low star formation rates lie consider- 
ably below L24 of normal giant elliptical galaxies. The 
black filled squares in Figure 8 illustrate the very small 
diffuse interstellar L24 emission for five representative el- 
lipticals computed with our steady-state dust evolution. 

6. CONCLUSIONS 

The banana-shaped infrared color-color correlation for 
normal galaxies - a plot oilog{Lio/ Lk) vs. log(_L24/iA') 

2 The two elliptical (red) galaxies with logL24 ~ 42.6 (NGC 
807 and UGC 01503) in Figure 8 have massive (~ 10^ Mq) multi- 
kpc rotating dusty HI disks (Goudfrooij et al. 1994; Young et 
al. 2009; Young 2002), likely to be old mergers. NGC 3125 (at 
logL24 ~ 42.26) is a compact dwarf galaxy (ct = 48 km s"-*^) with 
strong emission lines (Koprolin and Zeilinger 2000). Very little 
is known about NGC 4786 (logL24 ^ 42.26). The group of five 
E galaxies at logL24 ~ 41.92 includes NGC 5018 [with a large 
10* Mq filament of dusty HI gas], NGC 3557 [an FRI radio 
source (Birkinshaw & Davies 1985) with extended Ha emission 
(Goudfrooij et al. 1994)], NGC 4125 [with a central X-ray source 
(Flohic et al. 2006)], IC 4296 [an FRI radio galaxy with optical line 
emission and a dusty core], and IC 1459 [which contains detectable 
HI(Doyle et al. 2005)]. Finally, very little is known about the the 
outlier SO galaxy NGC 526 at logL24 » 43.45,logL70 ^ 43.0 in 
Figure 8. 
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(bottom panel of Figure 1) - approximately follows 
the morphological progression along Hubble types. As 
log(L24/L;f) and log(L7o/L_fs-) increase with the specific 
star-forming rate, the dominant Hubble type becomes 
progressively later: Sa,Sab Sb,Scd Irr. By con- 
trast, most E galaxies lie near the bottom of the ba- 
nana in a contiguous region of restricted \og{L24/ Lk) 
and decreasing \og{Ljo/LK)- Some gas-rich E galaxies 
lie along the star-forming locus of late type galaxies and 
SO galaxies are found throughout the banana. Our con- 
cern here is to understand the variation of \og{L2Al Lk) 
and \og{LjQ/ Lk) among normal elliptical galaxies in the 
banana plot. 

We show that the banana-shaped correlation, first dis- 
cussed by TBM09, becomes much tighter when galaxies 
with AGN emission are removed. Dust grains can be 
significantly heated by AGN radiation and emit strongly 
at 24/im within several arcseconds of the galactic cen- 
ters. These galaxies, having hard, highly concentrated 
infrared emission, can be identified by the presence of 
Airy diffraction rings at 24/im. AGN galaxies of all types 
with Airy rings lie above the infrared banana correlation. 
To achieve this, centrally illuminated dust would need 
to move along a iva.]eciOTy 5\og{L2i/ L k) / S\og{Lji^ / L k) 
that is steeper than the star formation trajectory along 
the upper banana, 5\og{L2i/ Lk) / 5\og{Ljo/ Lk) ~ 1-2. 
Because of this slope difference, we suggest that the con- 
centrated source of central IR radiation is due to AGN 
rather than small starburst regions unresolved at 24/im. 

An increasing specific star formation rate from early 
spiral to irregular galaxies with L24 oc iyo** apparently 
explains the progressive transition toward later galactic 
Hubble types upward along the banana. The bottom 
part of the infrared banana is populated by elliptical 
galaxies with a relatively small but significant scatter in 
24/im emission, \og{L24,/ Lk) ~ 30.15 ±0.15, and a much 
larger scatter in 70/im emission, 29 ^ log(L7o/i/<-) ^ 31 
in units of (erg s~^)/Lkq- 

We find that the scatter in \og{L24/ Lk) among el- 
liptical galaxies correlates with their stellar metallicity, 
{L24/LK) oc [Z/H]^-'^. Emission from hotter dust at 
24/im is likely to originate in circumstellar regions near 
AGB stars where dust is thought to be produced. The 
correlation of 24/im emission with metallicity is prob- 
ably due to the simple expectation that outfiows from 
metal-rich AGB stars produce more dust per gram of 
gas ejected. 

But the surface brightness variation in 24/im emis- 
sion from elliptical galaxies, while largely mimicking the 
de Vaucouleurs profile, contains unexpected color varia- 
tions when compared with the K-hand. When plotted 
against galactic radius normalized by the iiT-band effec- 
tive radius R/Re{K), the mean (K — 2A) color for a sam- 
ple of well-observed elliptical galaxies first decreases in 
a power-law fashion, then undergoes a broad minimum 
followed by a rise at large radius R/Re{K). The initial 
decrease is consistent with decreasing stellar metallicity 
with radius similar to that found at optical wavelengths, 
A[Z/H]/A\ogR « -0.3. The origin of the flattening 
and possible rise in (i^T — 24) beyond R/Re{K) ~ 1 is un- 
certain. If hot dust in AGB winds is thermally heated, in- 
creasing (K — 24) with decreasing stellar metallicity may 
be related to increasing effective temperature of AGB 



stars with decreasing metallicity. Alternatively, lower 
velocity (radiation pressure-driven) outflows from metal 
poor AGB stars might have larger 24/im emissivity. 

The intrinsic scatter of elliptical galaxies in the 
banana-shaped infrared correlation at longer wave- 
lengths, 70/im and 160/im, is much larger. Recent gas- 
rich mergers may be a source of cold dust for some ellip- 
tical galaxies. However, for most ellipticals we propose 
that this scatter is due to transient hydrodynamic flows 
of cold dust from galactic cores far out into the surround- 
ing hot galactic interstellar gas following AGN feedback 
events. The variation of log(L7o/ii<-) is uncorrelated 
with X-ray luminosity Lx of the galactic/group bary- 
onic halo. The Ljq/Lk ratios of some group-centered E 
galaxies like NGC 4472 and 1399 are in good agreement 
with a steady state model in which dust created by AGB 
stars is locally destroyed by sputtering at each radius in 
the galactic atmospheres. Other ellipticals such as NGC 
4636 and 5044 have much larger L70 luminosities from 
excess cold dust that cannot be produced by local stars. 
Lacking convincing evidence of recent gas-rich mergers 
in these otherwise normal elliptical galaxies, we propose 
that AGN outbursts cause excess dust to be hydrody- 
namically transported out to 5-10 kpc into the hot viri- 
alized gaseous atmospheres. Spatially extended 70/im 
emission has been confirmed in relatively nearby ellipti- 
cal galaxies (TBM07a,b). Evidently, values of Ljq/Lk 
for elliptical galaxies along the bottom of the banana 
plot are transient and may vary on ~ 10^ yr timescales, 
depending on the grain sputtering time and the AGN 
duty cycle. The outward relocation of variable masses of 
cold dust from the central reservoir, observed at different 
phases of the outflow, may explain why galaxies of sim- 
ilar Lk have L70 luminosities that vary over two orders 
of magnitude. 

However, the continuity of the banana plot from E 
galaxies to early spirals (where L24/LK and L-jq/Lk are 
due to star formation) requires that the enhanced 70/im 
emission from extended excess cold dust in E galaxies 
must not exceed L^q/Lk as produced by dust in star- 
forming galaxies of later types. To explain this unlikely 
coincidence, we speculate that star formation naturally 
begins in the cores of elliptical galaxies as the mass of 
accumulated cold gas and dust approaches some star- 
forming threshold. This conjecture follows from the 
expected secular accumulation of dust-cooled gas from 
AGB stellar mass loss in the central '^kpc of elliptical 
galaxies (Mathews & Brighenti 2003). 
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TABLE 1 

Galaxies in Figure 1 with Airy Diffraction'' 



Name 


T-Type'' 


D 


Log Lfcs 


LogL24/im 


LogZ/TO^m 






(Mpc) 




(erg s"^) 


(erg s~^) 


iNGC 0315 


-4 


58.9 


11.71 


42.68 


42.77 


M/"' no 1 /I 


-1.8 


20.7 


9.508 


42.65 


42.82 


iNCjO 05zo 


-2.0 


78.7 


■ ■ ■ 


43.44 


43.02 


InLtC izbo 


-2.1 


30.1 


10.51 


43.07 


43.77 


iNGG i3oD 


-0.8 


16.5 


10.56 


42.69 


42.98 


NGC 1377 


-2.1 


22.2 


10.15 


43.13 


43.2 


NGC 2110 


-3.0 


31.3 


11.08 


42.97 


43.4 


NGC 2325 


-4.7 


31.9 


11.2 


41.62 


41.3 


NGC 3226 


-4.8 


23.5 


10.66 


41.35 


41.97 


NGC 3265 


-4.8 


21.3 


9.83 


42.3 


42.8 


NGC 4261 


-4.8 


31.6 


11.44 


41.88 


41.81 


NGC 5077 


-4.8 


32.4 


11.08 


41.72 


41.85 


NGC 6776 


-4.1 


70.4 


11.43 


42.08 


42.5 


NGC 5273 


-1.9 


16.5 


10.32 


41.53 


41.96 


IC 5063 


-0.9 


45.3 


11.16 


43.82 


43.67 


E103-035 


-0.3 


54.2 


10.68 


43.87 


43.43 


ESO428-014 


-1.7 


21.3 


10.5 


42.99 


43.04 



^ These galaxies are identified with green squares in the central panel of Figure 1. 



Morphological type from HyperLeda. 
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TABLE 2 

Relevant Properties of the Tracer et al. Sample 



Name 


T-Type'' Log Lks 


LogL24^m 


LogLntim 


[Z/H]b 


[Z/H]b 


Age'' 


D R,{K) 




(NGC) 






(erg s~^) 


(■Re/8) 


(■Re/2) 


(Gyr) 


(Mpc) (//) 


(kpc) 



221 


-4.7 


9.13 


39.22 


37.93 


0.00±0.05 


-0.08± 0.05 


4.9±1.3 


0.81 


27.0 


0.11 


2.72 




-4.0 


11.71 


42.68 


42.77 


0.32±0.06 


0.17±0.06 


6.9±1.9 


58.88 


22.2 


6.36 


0.96 


507 


-3.3 


11.68 


42.02 




0.18±0.06 


0.28±0.13 


3.5±2.7 


67.19 


24.0 


7.84 


1.36 


584 


-4.6 


11.18 


41.48 


41.18 


0.48±0.03 


0.26±0.07 


3.4±1.1 


23.76 


23.3 


2.69 


1.62 


636 


-4.8 


10.67 


40.73 




0.34±0.07 


0.10±0.05 


6.8±1.4 


22.28 


15.3 


1.66 


1.39 


720 


-4.8 


11.14 


41.42 


40.69 


0.44±0.15 


1.13±0.42 


1.1±0.8 


22.29 


25.2 


2.73 


1.57 


821 


-4.8 


10.95 


41.06 




0.22±0.03 


0.12±0.05 


7.1±1.2 


24.09 


21.9 


2.56 


1.07 


1374 


-4.3 


10.68 


40.73 




0.13±0.07 




9.5±2.6 


19.77 


15.4 


1.48 


1.25 


1399 


-4.5 


11.40 


41.49 


40.50 


0.29±0.06 




11.5±2.4 


19.40 


32.9 


3.10 


2.83 


1404 


-4.7 


11.20 


41.45 


40.80 


0.25±0.05 




9.0±2.5 


19.40 


18.7 


1.76 


3.57 


1427 


-4.0 


10.66 


40.90 


40.89 


-0.07±0.03 




12.2±1.6 


23.55 


21.5 


2.46 


1.45 


1700 


-4.7 


11.27 


41.51 


41.35 


0.50±0.03 


0.32±0.05 


2.8±0.5 


38.04 


14.9 


2.76 


1.57 


2300 


-3.4 


11.18 


41.36 




0.36±0.04 


0.14±0.04 


10.1±2.0 


27.67 


22.3 


2.99 


1.54 


2778 


-4.7 


10.26 


40.37 




0.29±0.07 


-0.04±0.06 


14.9±3.5 


22.91 


10.5 


1.17 


1.67 


3377 


-4.8 


10.47 


40.50 


40.74 


0.19±0.06 


-0.12±0.04 


5.9±1.2 


11.22 


25.4 


1.39 


0.76 


3379 


-4.8 


10.89 


41.03 


40.56 


0.21±0.04 


0.00±0.04 


13.2±2.4 


10.57 


28.5 


1.46 


1.09 


3608 


-4.8 


10.83 


41.00 


41.27 


0.26±0.04 


0.07±0.07 


9.0±2.5 


22.91 


16.2 


1.80 


1.44 


4261'' 


-4.8 


11.44 


41.88 


41.81 


0.18±0.04 


-0.01±0.03 


21. Oil. 7 


31.62 


24.2 


3.72 


1.17 


4374 


-4.2 


11.39 


41.45 


42.03 


0.12±0.03 


-0.01±0.05 


14.4±2.8 


18.37 


33.5 


2.99 


1.12 


4472 


-4.8 


11.65 


41.71 


40.96 


0.25±0.05 


0.18±0.06 


8.4±2.7 


17.06 


56.1 


4.65 


0.89 


4478 


-4.8 


10.52 


40.74 




0.29±0.10 


-0.01±0.04 


10.3±1.7 


18.11 


11.8 


1.04 


1.35 


4489 


-4.8 


10.10 


40.11 




0.14±0.06 


-0.15±0.04 


4.6±0.5 


17.86 


14.1 


1.22 


1.13 


4552 


-4.6 


11.03 


41.26 


41.06 


0.27±0.04 


0.10±0.04 


13.1±2.7 


15.34 


22.8 


1.70 


1.65 


4649 


-4.6 


11.52 


41.58 


40.86 


0.27±0.04 


0.05±0.04 


18.3±2.8 


17.06 


42.1 


3.49 


0.46 


4697 


-4.8 


11.22 


41.26 


41.92 


0.06±0.05 


-0.29±0.04 


15.6±2.5 


16.22 


39.5 


3.11 


0.65 


5813 


-4.8 


11.40 


41.37 


41.51 


-0.04±0.03 


-0.20±0.06 


24.3±2.0 


32.21 


36.2 


5.66 


0.95 


5831 


-4.8 


10.84 


41.10 




0.54±0.04 


0.18±0.06 


4.2±1.0 


27.16 


19.4 


2.56 


0.75 


5846 


-4.7 


11.36 


41.43 


41.53 


0.14±0.05 


-0.13±0.04 


23.1±2.7 


24.89 


32.7 


3.96 


2.47 


6703 


-2.8 


11.06 


41.34 


41.36 


0.30±0.06 


0.03±0.06 


7.9±2.1 


32.06 


20.1 


3.13 


1.16 


7619 


-4.7 


11.58 


41.62 




0.20±0.03 


0.07±0.04 


15.1±2.2 


52.97 


16.4 


4.22 


1.64 


7626 


-4.8 


11.34 


41.44 


41.17 


0.16±0.04 


-0.06±0.02 


18.2±2.0 


39.99 


20.7 


4.02 


1.13 


7785 


-4.8 


11.46 


41.58 


41.45 


0.20±0.01 


-0.01±0.04 


15.0±2.3 


55. 16 


17.2 


4.64 


0.84 



^ Morphological type from HyperLeda. 

^ Values of [Z/H] for Re/8 and Re/2 are "model 4" values from Trager et al. (2000a, b) refer to the B-band effective radius Re which on average is about 1.45 times 
larger than the K"-band value Re(i^) (TBM08). Ages listed above refer to the -Re/2 observations. 
^ Maximum extent of observable 24fj,ni emission. 

Identified as an AGN galaxy from 24/j,m Airy rings. 
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Fig. 1. — Upper panel: Infrared color-color banana plot L24/LK vs. L-jq/Lk for SINGS galaxies showing progression of Hubble types. 
Central panel: IR color-color plot showing galaxies in the TBM09b sample with X indicating SINGS galaxies. Green squares indicate 
TBM09b (AGN or starburst) galaxies with centrally concentrated 24/xm emission. Bottom panel: Same as central panel without galaxies 
having centrally concentrated 2ifim emission. 
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Fig. 2. — Upper panels: 24/jm images of NGC 1398 and 315 showing Airy diffraction rings and NGC 4472 with diffraction-free image. 
The colors and contours are qualitative only, intended to illustrate non-axisymmetric features in the diffraction pattern. Lower panel: 
Azimuthally averaged radial diffraction fluxes at 24/^m for the three galaxies imaged above: NGC 1386 (black dotted line); NGC 0315 
(blue dashed line); NGC 4472 (light blue long dash-triple dotted line). For comparison we plot the point response function of the Spitzer 
telescope at 24/jm (green solid line). All fluxes are normalized to 100 (in arbitrary units) at one arcsecond. 
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Fig. 3. — Expanded view of E and E-SO galaxies (with T < —2.7) at the bottom of the banana plot, omitting AGN galaxies with 
significant point source emission at 24^m. Red circles identify a subset of these galaxies that are also in the sample of Trager et al. (2000a, 
b) which have been observed at both 24 and 70/im but omitting two that are identified as AGNs in Table 1, NGC 0315 and NGC 4261. 
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Fig. 4. — Correlation between L2a/Lk and stellar metallicity [Z/H\ for the Trager et al. (2000a, b) sample of E galaxies. Upper and 
lower panels show the linear relation with metallicities values derived in the Re/S and Re/2 aperture respectively where Re is the effective 
radius in the B-band. Two galaxies identified as AGNs in Table 1 (NGC 315 and NGC 4261) are not included in either panel. In the lower 
panel we also exclude the outlying (and very uncertain) value [Z/H] = 1.13 for NGC 720. 
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Fig. 5. — Upper panel: Variation of Trager sample-averaged radial color profile (in surface brightness units) with galaxy radius normalized 
to the K-band half light radius Re{K). Lower panel; Variation of Trager sample-averaged radial stellar metallicity profiles found using the 
correlation in Figure 4 derived for the ile/2 aperture. 
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Fig. 6. — Plot of L70 vs. Lk for elliptical galaxies (T < —3.5) in the TMB09b sample. Galaxies marked in red all have excess extended 
dust and the large red circles indicate galaxies in which the spatial extension has been observed at 70/im. Vertical green arrows show the 
approximate estimated increase in L70 due to excess cold dust in NGC 5044 and 4636. 
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Fig. 7. — Infrared color-color plot as in Figure 3 showing only E and E-SO galaxies with T < —3.5. Galaxies in red and circled in red are 
identical with those in Figure 6, but now appear along the right boundary of the lower banana. As in Figure 6 the horizontal green arrows 
show the trajectories of NGC 5044 and 4636 when extended cold dust is added to 70/im emission produced only by local AGB stellar mass 
loss. 
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Fig. 8. — Plot of L24 against L70 as in the bottom panel of Figure 1 but without the Lk normalization. Symbols are identical to 
those in that panel. Filled squares show computed steady-state emission from diffusely distributed interstellar dust in five representative 
elliptical galaxies (NGC 4472, 1404 4636, 1399 & 5044 in order of decreasing L24)- The star formation locus of SINGS galaxies (crosses) is 
L24, oc Li^OS. 



